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Beschreibung 

[0001] Die Erfindung betrifft ein Verfahren zum modellgestutzten Bestimmen der in die Zylinder einer Brennkraftma- 

schine einstromenden Luftmasse nach dam Oberbegriff des Patentanspruches 1 . 

5 [0002] Motorsteuerungssysteme fur Brennkraftmaschinen, die mit Kraftstoffeinspritzung arbeiten, benotigen die vom 
Motor angesaugte Luftmasse m2y/ als ein MaB fur die Motorlast. Diese Gr6f3e bildet die Basis zur Realisierung eines 
geforderten Kraftstoff-Luft-Verhaltnisses. Wachsende Anforderungen an Motorsteuerungssysteme, wie die Verringe- 
rung der Schadstoffemission von Kraftfahrzeugen, bedingen, da3 die LastgroBe fur stationare und instationare Vor- 
gange mit geringen zulassigen Fehlern bestimmt werden mu3. Neben den gehannten Betriebsfallen bietet die genaue 

10 Lasterfassung wahrend der Warmlaufphase der Brennkraftmaschine ein erhebliches Potential zur Schadstoffreduktion. 
[0003] Bei luftmassengefuhrten Motorsteuerungssystemen stellt im Instationarbetrieb das als Lastsignal der Brenn- 
kraftmaschine dienende Signal des Luftmassenmessers, der stromaufwarts des Saugrohrs angeordnet ist, kein Ma3 
fur die tatsachliche Fullung der Zylinder dar, weil das Volumen des Saugrohrs stromabwarts der Drosselklappe als 
Luftspeicher wirkt, der befullt und entleert werden muB. Die maBgebende Luftmasse fur die Einspritzzeitberechnung 

IS ist aber diejenlge Luftmasse, die aus dem Saugrohr heraus und in den jeweiligen Zylinder hineinstromt. 

[0004] Bei saugrohrdruckgefuhrten Motorsteuerungssystemen gibt zwar das Ausgangssignal des Drucksensors die 
tatsachlichen Druckverhaltnisse im Saugrohr wieder. die MeBgroBen stehen aber u.a. aufgrund der notwendigen Mit- 
telung der MeBgroBe erst reiativ spat zur Verfugung. 

[0005] Mit der Einfuhrung variabler Ansaugsysteme und variabler Ventilsteuerungen entstehen fur empirisch gewon- 
20 nene Modelle zur Gewinnung der LastgroBe aus MeBsignalen eine sehr groBe Vielzahl von EinfluBgroBen, die die 
entsprechenden Modellparameter beeinflussen. 

[0006] Auf physikalischen Ansatzen basierende modellgestutzte Berechnungsmethoden stellen einen guten Aus- 
gangspunkt zur genauen Bestimmung der Luftmasse m^y/ dar. 

[0007] Aus der DE 39 1 9 488 C2 ist eine Vorrichtung zur Regelung und zur Vorausbestimmung der Ansaugluftmenge 
25 einer saugrohrdruckgefuhrten Brennkraftmaschine bekannt, bei der der Drosselklappenoffnungsgrad und die Motor- 
. drehzahl als Grundlage zur Berechnung des derzeitigen Wertes der in den Brennraum der Maschine eingesaugten 
Luft verwendet werden. Diese berechnete, gegenwartige Ansaugluftmenge wird dann als Grundlage zur Berechnung 
des vorausbestimmten Wertes fur die Ansaugluftmenge, die in den Brennraum der Maschine zu einer bestimmten Zeit 
von dem Punkt an, an dem die Berechnung ausgefOhrt wurde, einzusaugen ist, benutzt. Das Drucksignal, das strom- 
30 abwarts der DrosselklajDpe gemessen wird, wird mit Hilfe von theoretischen Beziehungen korrigiert, so daB eine Ver- 
besserung der Bestimmung der angesaugten Luftmasse erreicht und damit eine genauere Berechnung der Einspritz- 
zeit moglich ist. 

[0008] Im instationaren Betrieb der Brennkraftmaschine ist es aber wunschenswert, die Bestimmung der in die Zy- 
linder einstromenden Luftmasse noch genauer durchzufuhren. 
35 [0009] Der Erfindung liegt die Aufgabe zugrunde ein Verfahren anzugeben, mit dem die tatsachlich in den Zylinder 
der Brennkraftmaschine einstromende Luftmasse mit hoher Genauigkeit bestimmt werden kann. AuBerdem sollen 
systembedingte Totzeilen, die aufgrund der Krafts toffvorlagerung und der Rechenzeit bei der Berechnung der Ein- 
spritzzeit auftreten konnen, kompenslert werden. 

[0010] Diese Aufgabe wird gemaB den Merkmalen des Patentanspruches 1 gelost. 

40 [0011] Vorteilhafte Weiterbildungen finden sich in den Unteranspruchen. 

[0012] Ausgehend von einem bekannten Ansatz ergibt sich eine Modellbeschreibung, die auf einer nichtlinearen 
Differentialgleichung basiert. Im folgenden wird eine Approximation dieser nichtlinearen Gleichung vorgestellL Im Er- 
gebnis dieser Approximation laBt sich das Systemverhalten mittels einer bilinearen Gleichung beschreiben, die die 
schnelle Losung der Beziehung im Motorsteuergerat des Kraftfahrzeugs unter Echtzeitbedingungen gestattet. Der 

45 gewahlte Modellansatz beinhaltet dabei die Modellierung von variablen Saugsystemen und Systemen mit variablen 
Ventilsteuerungen. Die durch diese Anordnung und durch dynamische Nachladung, d.h. durch Reflexionen von Druck- 
wellen im Saugrohr hervorgerufenen Effekte, konnen ausschlieBlich durch die Wahl stationar bestimmbarer Parameter 
des Modelles sehr gut berOcksichtigt werden. Alle Modellparameter sind einerseits physikalisch interpretierbar und 
andererseils ausschlieBlich aus stationaren Messungen zu gewinnen. 

so [0013] Die melslen Algorithmen zur zeildiskreten Losung der Differentialgleichung, die das Verhallen des hier ce- 
nutzten Modelles beschreibt, erfordern vor allem bei geringem Druckabfail Ober der Drosselklappe, d.h. bei Vollast 
eine sehr kleine Rechenschrillweite, um numerisch stabil zu arbeiten. Die Folge ware ein unvertretbarer Rechenaut- 
wand bei der Bestimmung der LastgroBe. Da Lasterfassungssysteme meist segmentsynchron arbeiten, d.h. fur 4-Zy- 
lindermotoren wird alle 180° KW ein MeBwert abgetastet, muB die Modellgleichung ebenfalls segmentsynchron gelost 

55 werden. Im nachfolgenden wird ein absolut stabiles Differenzenschema zur Losung von Differentialgleichungen ein- 
gesetzt, das numerische Stabilitat bei beliebiger Schrittweite garantiert. 

[001 4] Das erfindungsgemaBe modellgestutzte Berechnungsverfahren bietet zudem die Moglichkeil einer Pradiklion 
des Lastsignales um eine wahlbare Anzahl von Abtastschritten, d.h. eine Vorhersage des Lastsignales mit variablem 
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Pradlkliohshorizont. Wird die dem Pradiktionshorizont bei konstanter Drehzahl proporlionale Pradiktionszelt nicht zu 
gro3, so erhall man ein pradiziertes Lastsignal hoher Genauigkeit. 

[001 5] Eine solche Vorhersage ist notwendig, da zwischen der Erfassung relevanter MeBwerte und der Berechnung 
der LastgroBe eine Totzeit entsteht. Desweiteren muB aus Grunden der Gemischaufbereitung vor dem eigentlichen 

s Beginn der Ansaugphase des jeweiligen Zytinders moglichst genau die Kraftstoffmasse uber die Einspritzventite zu- 
gemessen werden. die im Verlauf der kommenden Ansaugphase im gewunschten Verhaltnis zur Luftmasse m^^/Steht. 
Ein variabler Pradiktionshorizont verbessert die Gute der Kraftstoffzumessung im instationaren Motorbetrieb. Da bei 
steigender Drehzahl die Segmentzeit abnimmt, muB der Einspritzvorgang eine groBere Anzahl von Segmenten eher 
beginnen, als dies bei einer niedrigeren Drehzahl der Fall ist. Um die zu dosierende Kraftstoffmasse moglichst exakt 

10 bestimmen zu konnen, ist die Pradiktion der LastgroBe um die Anzahl von Segmenten, um die die Kraftstotfvorlagerung 
vorgenommen wird, notwendig, um ein gefordertes Kraftstoff-Luft -Verhaltnis auch in diesem Fall einzuhalten. Die Pra- 
diktion der LastgroBe tragt somit aus einer wesentlichen Verbesserung der Einhallung des gefordertene Kraftstoff- 
Luft-Verhaltnisses im instationaren Motorbetrieb bei. Dieses System zur modellgestutzten Lasterfassung ist. in den 
bekannten Motorsteuerungssystemen, d.h. bei luftmassengefuhrte bzw. saugrohrdruckgefuhrte Motorsteuerungssy- 

IS Sterne wird im folgenden ein Korrekturalgorithmus in Form eines Modellregelkreises formuliert, der bei auftretenden 
Ungenauigkeiten von Modellparametern eine permanente Genauigkeitsverbesserung, d.h. einen Modellabgleich im 
stationaren und instationaren Betrieb gestattet. 

[0016] Ein Ausfuhrungsbeispiel des erfindungsgemaBen Verfahrens wird anhand der nachfolgenden schematischeh . 
Zeichnungen beschrieben. Dabei zeigen: 

20 

Figur 1 eine Prinzipskizze zum Saugsystem einer Otto-Brennkraftmaschine einschlieBlich der entsprechenden Mo- 
dell- und MeBgroBen, 
Figur 2 die DurchfluBfunktion und die dazugehorige Polygonzugapproximation, 

Figur 3 eine Prinzipdarstellung zum Modellregelkreis fur luftmassengefuhrte Motorsteuerungssysteme und 
25 Figur 4 eine Prinzipdarstellung zum Modellregelkreis fur saugrohrdruckgefuhrte Motorsteuerungssysteme. 

[0017] Bei der modellgestutzten Berechnijng der LastgroBe /TJ^y/ wird von der in Figur 1 dargestellten prinzipiellen 
Anordnung ausgegangen. Aus Grunden der Ubersichllichkeit wird dabei nur ein Zylinder der Brennkraftmaschine dar- 
gestellt. Mit dem Bezugzeichen 10 Ist dabei ein Saugrohr einer Brennkraftmaschine bezeichnet, in. dem eine Drossel- 

30 klappe 11 angeordnet ist. Die Drosselklappe 11 ist mit einem, den Offnungsgrad der Drosselklappe ermittelnden Dros- 
selklappenstellungsfuhler 14 verbunden. Stromaufwarts der Drosselklappe 11 ist bei einem luftmassengefuhrten Mo- 
torsteuerungssystem ein Luftmassenmesser 12 angeordnet, wahrend bei einem saugrohrdruckgefuhrten Motorsteue- 
rungssystem ein Saugrohrdruckfuhler 13 im Saugrohr angeordnet 1st. Je nach Art der Lasterfassung ist somit nur eine 
derbeiden Komponenten 12, 13vorhanden. Die Ausgange des Luftmassenmessers 12, des Drosselklappenstellungs- 

35 gebers 14 und des zum Luftmassenmesser 12 alternativ vorhandenen Saugrohrdrucksensors 13 sind mit Eingangen 
einer nicht dargestellten , an sich bekannten elektronischen Steuerungseinrichtung der Brennkraftmaschine verbunden. 
AuBerdem sind in Figur 1 noch schematisch ein EinlaBventil 15, ein AuslaBventil 16, sowie ein in einem Zylinder 17 
beweglichen Kolben 18 dargestellt. 

[0018] AuBerdem sind in Figur 1 ausgewahlte GroBen bzw. Parameter des Saugsystems eingezeichnet. Dabei be- 
40 deutet das Dachsymbol uber einer GroBe, daB es sich um eine ModellgroBe handelt, wahrend GroBen ohne Dach- 
syrhboi MeBgroBen reprasentieren. Im einzelnen bedeuten: 

[0019] Py Umgebungsdruck. Ps Saugrohrdruck, Ts Temperatur der Luft im Saugrohr, Vg das Volumen des Saug- 
rohrs. 

[0020] GroBen mit einem Punktsymbol kennzeichnen die erste zeitliche Ableitung der entsprechenden GroBen. 
45 isi somit der Luftmassenstrom an der Drosselklappe und m^yi ist der Luftmassenstrom der tatsachlich in den Zylinder 

der Brennkraftmaschine einstromt. 

[0021] Die grundlegende Aufgabe bei der modellgestutzten Berechnung des Motorlastzusta'ndes besteht nun in der 
Losung der Differentialgleichung fur den Saugrohrdruck 



so 



55 



J 



(2.1) 



die sich unter der Voraussetzung konstanter Temperatur der Luft im Saugrohr Tg aus der Zustandsgleichung Idealer 
Gase herleiten Ia3t. 

[0022] Mit Rl ist dabei die allgemeine Gaskonstante bezeichnet. 
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[0023] Die LastgroBe m2y/ wird durch Integration aus denn Zylindernnassenstromiii^ri bestimmt. Die durch (2.1) be- 
schriebenen Verhaltnisse sind auf Mehrzylinder-Brennkraftmaschinen mit Schwingrohr- (Schaltsaugrohr-) und/oder 
Resonanzsaugsysteme ohne strukturelle Anderungen anwendbar. 

[0024] Fur Systeme mit Multi-Point-Einspritzungen, bei denen die Kraftstoffzumessung durch mehrere Einspritzven- 
tile erfolgt, gibt die Gleichung (2.1) die Verhaltnisse genauer wieder als dies bei Single-Point-Einspritzungen, d.h. bei 
Einspritzungen, bei denen der Kraftstoff mittels eines einzigen Kraftstoffelnspritzventiles zugemessen wird, der Fall 
ist. Bei erstgenannter Art der Kraftstoffzumessung ist nahezu das gesamte Ansaugsystem mit Luft gefullt. Lediglich in 
einem kleinen Bereich vor den EinlaBventilen befindet sich ein Kraftstoff-Luftgemisch. Im Gegensatz dazu ist bei Single- 
Point-Einspritzsystemen das gesamte Saugrohr von der Drosselklappe bis zum EinlaBventil mit Kraftstoff-Luft-Ge- 
misch gefullt, da das Einspritzventil vor der Drosselklappe angeordnet Ist. In diesem Fall stellt die Annahme eines 
idealen Gases eine starkere Naherung dar, als dies bei der Multi-Point-Einspritzung der Fall ist. Bel Single-Point- 
Einspritzung erfolgt die Kraftstoffzumessung entsprechendmo* bei Multi-Point-Elnspritzung entsprechend^^^ 
[0025] Im folgenden wird die Berechnung der Massenstrome moK und m^w naher beschrieben. 
[0026] Die ModellgroBe des Luftmassenstromes an der Drosselklappe mm- wird durch die DurchfluBglelchung idealer 
Gase durch Drosselstellen beschrieben. An der Drosselstelle auftretende Stromungsverluste werden durch den redu- 
zierten Stromungsquerschnitt Apf^QberucksichUgX. Der Luftmassenstrom wird demnach durch die Beziehung 



20 
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fur uberkritische Druckverhaltnisse bzw. 



y = const, fur kritische Druckverhaltnisse 



(2.2) 



bestimrrit. 



45 



so 



moK ■ 

k: 
Rl: 

Is: 



ModellgroBe des Luftmassenstromes an der Drosselklappe 

reduzierter Stromungsquerschnitt 

Adiabatenexponent 

allgemeine Gaskonstante 

Temperatur der Luft im Saugrohr 

ModellgroBe des Umgebungsdruckes 

ModellgroBe des Saugrohrdruckes 

DurchfluBfunktion. 



[0027] An der Drosselstelle, d.h. an der Drosselklappe auftretende Stromungsverluste werden uber die geeignete 
Wahl von A^^q berucksichtigt. Aus stationaren Messungen kann bei bekannten Drucken vor und hinter der Drosselst- 
elle und bekanntem Massenstrom durch die Drosselstelle eine Zuordnung zwischen dem vom Drosselklappenstel- 
55 lungsfuhler 14 ermittelten Drosselklappenwinkel und dem entsprechendem reduzierten Querschhitt A^^^ angegeben 
werden. 

[0028] Wird der Luftmassenstrom an der Drosselklappe durch die Beziehung (2.2) beschrieben. so entsteht ein 
komplizierter Algorithmus zur numerisch richtigen Losung der Differentialgleichung (2.1). Zur Reduktion des Rechen- 



45 
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aufwandes wird die DurchfluBfunktion y durch einen Polygonzug approximiert. 

[0029] FIgur 2 zeigt den Verlauf der DurchfluBfunktion >|/ und das darauf angewandte Approximationsprinzip. Inner- 

halb eines Abschnittes i (1 = 1...k) wird die DurchfluBfuntion y durch eine Gerade dargestellt. Mit einer vertretbaren 
Anzahl von Geradenabschnitten kann damit eine gule Approximation erreicht werden. Durch einen solchen Ansatz 
kann die Gleichung (2.2) zur Berechnung des Massenstromes an der Drosselklappe durch die Beziehung 



10 mOK APPROX = ^RED -\ ^ \ ' ' ' Pu 



I 2k 1 " 



(2.3) 



furi = (l...k) 
IS approximiert werden. 

[0030] In dieser Form beschreibt mj die Steigung und nj das Absolutglied des jeweiligen Geradenabschnittes. Die 
Werte fur die Steigung und fur das Absolutglied werden in Tabellen als Funktion des Verhaltnisses Saugrohrdruck zu 
Umgebungsdruck-^ abgelegt.' ' p 

[0031] Auf der, Ab^zisse von Figur 2 ist dabei das Druckverhaltnis -^^ und auf der Ordinate der Funktiohswert (0 - 
20 0.3) der DurchfluBfunktion \\f aufgetragen. ^ 
[0032] Fur Druckverhaltnlsse 

! 

25 P^'VK^^) 

ist Y = konstant, d.h. daB der DurchfluB an der Drosselstelle nur noch vom Querschnitt abhangig ist und nicht mehr 
von den Druckverhaltnissen. Die in die jeweiligen Zylinder der Brennkraftmaschine einstromende Luftmasse laBt sich 
analytisch nur schwer bestimmen, da sie stark vom Ladungswechsel abhangt. Die Fullung der Zylinder wird weitest- 

30 gehend durch den Saugrohrdruck, die Drehzahl und durch die Ventilsteuerzeiten bestimmt. 

[0033] Zur moglichst genauen Berechnung des Massenstroms in den jeweiligen Zylinder «zy. ist deshalb einerseits 
die Beschreibung der Verhaltnisse im Ansaugtrakt der Brennkraftmaschine mittels partieller Differentialgleichungen 
und andererseits die Berechnung des Massenstromes am EinlaBventil nach der DurchfluBgleichung als erforderliche 
Randbedingung notwendig. Erst dieser komplizierte Ansatz gestattet die Berucksichtigung dynamischer Nachladeef- 

35 fekte, die von der Drehzahl, der Saugrohrgeometrie, der Zylinderzahl sowie den Ventilsteuerzeiten maBgeblich beein- 
fluBt werden. 

[0034] Da eine Berechnung nach oben genanntem Ansatz in der elektronischen Steuerungseinrichtung der Brenn- 
kraftmaschine n^cht realisierbar ist, geht eine mogliche Naherung von einem einfachen Zusammenhang zwischen 
Saugrohrdruck Pq und Zyindermassenstrom l^., aus. Fur einen weiten Bereich der sinnvollen Ventilsteuerzeiten kann 
40 dafur in guter Naherung von einem linearen Ansatz der Form' 



mzy/_APPROX =yi ' Ps-^Jo (2.4) 



ausgegangen werden. 

[0035] Die Steigung Yi und das Absolutglied Yq der Beziehung (2.4) sind dabei, unter Berucksichtigung aller wesent- 
lichen EinfluBfaktoren Funktionen der Drehzahl, der Saugrohrgeometrie, der Zylinderzahl, der Ventilsteuerzeiten sowie 

50 der Temperatur der Luft im Saugrohr Tg. Die Abhangigkeit der Werte von und Yq von den EinfluBgroBen" Drehzahl, 
Saugrohrgeometrie, Zylinderzahl und den Ventilsteuerzeiten und Ventilerhebungskurven kann dabei uber slationare 
Messungen ermittelt werden. Uber diese Wertebestimmung wird ebenfalls der EinluB von Schwingrohr- und/oder Re- 
sonanzsaugsystemen auf die von der Brennkraftmaschine angesaugte Luftmasse gut wiedergegeben. Die Werte von 
Y-i und Yo sind in Kennfeldern der elektronischen Motorsteuerungseinrichtung abgelegt. 

55 [0036] Als bestimmende GroBe zur Ermittlung der Motorlast wird der Saugrohrdruck ausgewahlt. Mil Hilfe der 
Modell-Differentialgleichung-soll diese GroBe moglichst exakt und schnell geschatzl werden. Die Schalzung von Pq 
erfordert die Losung der Gleichung (2. l ). 

[0037] Mit den anhand der Formein (2.2) und (2.3) eingefuhrten Vereinfachungen kann (2.1 ) durch die Beziehurig 
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(2.5) 



10 



IS 



approximiert werden. Betrachtet man, entsprechend den Voraussetzungen zur Herleitung von Gleichung (2.1), die 
Temperatur der Luft im Saugrohr Tg als eine langsam veranderliche MeBgroBe sowis Af^^Q als Eingangsgr63e, so 
laBt sich die nichtlineare Form der Differentialgleichung (2.1) durch die bilineare Gleichung (2,5) approximleren. 
[0038] Zur Losung der Gleichung (2.5) wird diese Beziehung in eine geeignete Differenzengleichung ubergefuhrt. 
[0039] Als Kriterium zur Auswahl des geeigneten Differenzenschemas konnen die folgenden prinzipiellen Anforde- 
rungen an die Losungseigenschaften der zur bildenden Differenzengleichung formuliert werden: 



20 



1. Das Differenzenschema muB auch unter extremen dynamischen Anforderungen konservativ sein, d.h. die Lo- 
sung der Differenzenglelchnung muB der Losung der Differentialgleichung entsprechen, 

2. die numerische Stabilitat muB zu Abtastzeiten, die den maximal moglichen Segmentzeiten entsprechen, im 
gesamten Arbeitsbereich das Saugrohrdruckes garantiert sein. 



[0040] Forderung 1 ist durch einen impliziten Rechenalgorithmus erfullbar. Aufgrund der Approximation der nichtli- 
nearen Differentialgleichung (2.1) durch eine bilineare Gleichung ist das entstehende implizite Losungsschema ohne 

25 ■ Einsatz iterativer Verfahren Ipsbar da die Differenzengleichung in eine explizite Form uberfuhrt werden kann. 

[0041] Die zweite Forderung ist aufgrund der Konditionierung der Differentialgleichung (2.1) und deren Approxima- 
tion (2.5) nur durch eine Rechenvorschrift zur Blldung der Differenzengleichung erfullbar. die absolut stabil arbeitet. 
Diese Verfahren werden auch als A-stabile Verfahren bezeichneL Kennzeichnend fur diese A-Stabilitat ist die Eigen- 
schaft des Algorithmus, bei einem stabilen Ausgangsproblem fur beliebige Werte der Abtastzeit, d.h. Segmentzeit T^j^ 

30 numerisch stabil zu sein. Eine mogliche Rechenvorschrift zur numerischen Losung von Differentialgleichungen, die 
beiden Forderungen gerecht wird, ist die Trapezregel. 

[0042] Die durch Anwendung der Trapezregel entstehende Differenzengleichung lautet im vorliegenden Fall 



35 



40 



fiir N = (1 . , CC7] 



(2.6) 



definiert. 

[0043] Wird diese Vorschrift auf (2.5) angewandt, so ergibt sich die Beziehung 
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fur N = {l...<5t?) und i = (l...k) 



(2.7) 



zu'r Berechnung des Saugrohrdruckes PslN] als Ma3 fur die Motorlast. 

[0044] [N] bedeutet dabei das aktuelle Segment bzw. der aktuelle Rechenschritt, [N+^ ].das nachstfolgende Segment 
bzw. dernachstfolgende Rechenschritt. 

[0045] Irri folgenden wird die Berechnung des^aktueKen und pradizierten Lastsignales beschrieben. 
[0046] Aus dem berechneten Saugrohrdruck P^kann der Luftmassenstrom n,2y, der in die Zylinder einstromt, durch 
die Beziehung (2.4) ermittelt werden. Wendet man einen einfachen Inte.grationsalgorithmus an, so erhait man fur die 
wahrend eines Ansaugtaktes von der Brennkraftmaschine angesaugte Luftmasse die Beziehung * - 
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mzy! [A^-l] + 'WZv/ [N] 
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fur N = (1 . . jCP) 



(2 . 8) 



40 



45 



[0047] Dabei wird davon au.sgegangen, daG der Anfangswert der LastgroBe null ist! Fur die' segmentsynchrone La- 
sterfassung sink! mit steigender Drehzahl die Segnnentzeit, wahrend die Segmentanzahl, um die eine Kraftstoffvorla- 
cei-vng vorgenommen wird, steicen rnuB. Aus diesem'Grund ist es erforderlich, die Pradiktion des Lastsignals fur ernen 

veranrieriichen Pradiktionshorizont H, d.h. fur eine bestimmte, in ersteV Linie drehzahlabhangige Anzahl H von Seg- 
rnerdhCi, auszulegen. Berucksichiigl rnsn diesen veranderiichen Piadikiionshorizonl H, so kann Gleichung (2.8) in der 
Fom': 
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mzyl 

V 

fiir N ^ ( 1 . . ) 



(2 . 9) 
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geschrieben werden. 

[0048] Fur die weiteren Uberlegungen wird davon ausgegangen, daG sich die Segmentzeit und die Parameter 
7-j und Yo der Beziehung (2.4); die zur Bestimmung des Massenstromes I/,^.aus dem Saugrohrdruck erforderlich 
sind, Qber die Pradiktionszeit nicht andern. Unter dieser^Voraussetzung wird die Pradiktion eines Wertes fur mz,,[N+H\ 
durch die Pradiktion des enlsprechenden Druckwerles PslN+H] erreicht. Dadurch nimmtdie Gteichung (2.9) die Form 
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mzyi\ 

fur N = (1 . . .CP) (2 . 10) 



an. 

A 

[0049] Da bei dem beschrlebenen Verfahren die zeitliche Anderung des Saugrohrdruckes Pq in analytischer Form 
vorliegt, wird im folgenden die Pradiktion des Druckwertes Pq [N+H\ durch H-fache Anwendung der Trapezregel er- 
reicht. In diesem Fall erhalt man die Beziehung 



fur N = (1 . . PC> ) (2 . 11) . 

A 

[0050] Bestimmt man den Druck [A/+/-/-1 ] in analoger Weise, so kann fur das pradizierte Lastsignai die Gleichung 



fur N = (1 . .00) (2 .12) 



angegeben warden. 

[0051] Wahit man fur den Pradiktionshorlzont H Werte In der GroIBen-ordnung von 1 ...3 Segmenten, so kann mit der 
Formel (2.12) ein gut pradiziertes Lastsignai erhalten werden. 

[0052] Im folgenden wird das Prinzip des Modellabgleichs tur luttmassen- und saugrohrdruckgetuhrte Motorsteue- 

rungssysteme erklart. 

[0053] Bedingt durch den EInsatz von Motoren mit variabler Ventilsteuerung und/oder veranderlicher Saugrohrgeo- 
metrie, durch Fertigungslolefanzen und Alterungserscheinungen, sowie durch Temperatureinflusse sind die Werte von 
Yi und Yo mit einer gewissen Unsicherheit behaftet. Die Parameter der Gleichung zur Bestimmung des Massenstromes 
in den Zylindern sind, wie oben beschrieben, Funktionen vielfaltiger EinfluBgroBen, von denen nur die wichtigsten 
erfaBt werden konnen. 

[0054] Bei der Berechnung des Massenstromes an der Drosselklappe wirken sich MeBfehler bei der Erfassung des 
Drosselklappenwinkels und Approximationsfehler bei der Polygonzugapproximation der Durchflu3funktion \\f auf die 
Modeilgro3en aus. Besonders bei kleinen Drosselklappenwinkein Ist die Systemempfindlichkeit gegenuber erstge- 
nannten Fehlern besonders hoch. Daraus ergibt sich, daB kleine Anderungen der Drosselklappenstellung einen gra- 
vierenden EInluB auf Massenstrom bzw. Saugrohrdruck haben. Um die Wirkung dieser Einflusse zu reduzieren, wird 
im folgenden ein Verfahren vorgeschlagen, das es gestattet, bestimmte GroBen, die EinfluB auf die Modellrechnung 
haben, so zu korrigieren, daB eine genaulgkeitsverbessernde Modellanpassung fur stationaren und instationaren Mo- 
torbelrieb durchgefuhrt werden kann. 

[0055] Die Anpassung wesentlicher Parameter des Modells zur Bestimmung der LastgroBe der Brennkraftmaschine 
erfolgt durch die Korrektur des aus dem gemessenen Drosselklappenwinkel bestimmten reduzierten Querschnltts A^j^q 
durch die KorrekturgroBe AAp^^p. 

[0056] Die EingangsgroBe zur korrigierten Saugrohrdruckberechnung Apj^Q wird damll durch die Beziehung 
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^REDKORR'^RED ^^RED (311) 



beschrieben. 

[0057] In der Gleichung (2.2) und nachfolgenden Formein wird dann A^^^ durch ^redkorr ©fsetzt. Zur Verbesse- 
rung des Folgeverhaltens des Regefkreises wird der aus dem MeBwert des Drosselklappenwinkels abgeleitete redu- 
zierte Drosselklappenquerschnitt A^^Qindie Modellrechnung einbezogen. Die KorrekturgroBe AA^ry^j wird durch Rea- 
lisierung eines Modellregelkreises gebildet. 

[0058] Fur luftmassengef uhrte Motorsteuerungssysteme ist der mittels des Luftmassenmessers an der Drosselklap- 
pe gernessene Luftmassenstrom rnoK^LMM FuhrungsgroBe dieses Regelkreises, wahrend fur saugrohrdruckge- 
f uhrte Systeme der gernessene Saugrohrdruck P^als FOhrungsgroBe genutzt wird. Uber eine Folgeregelung wird der 
Wert von AA/^^qSo bestimnnt, daB die Regelabweichung zwischen FOhrungsgroBe und der ensprechenden RegelgroBe 
minimiert wird. 

[0059] Urn auch im dynannischen Betrieb Genauigkeitsverbesserungen mit der genannten Methode zu erreichen, 
muB die MeBwerterfassung der FOhrungsgroBe nnoglichst exakt nachgebildet werden. In den meisten Fallen sind dabei 
das dynamische Verhalten des Sensors, d.h. entweder des Luftmassenmessers oder des Saugrohrdrucksensors und 
eine nachfolgend durchgefuhrte Mittelwertbildung zu berucksichtigen. 

[0060] Das dynamische Verhalten des jeweiligen Sensors kann in erster Naherung als ein System erster Ordnung 
mit eventuell arbeitspunktabhangigen Verzogerungszeiten modelliert werden. Im Falle eines luftmassengefOhrten 
Systems lautet eine mogliche Gleichung zur Beschreibung des Sensorverhaltens 
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fTlDK LMM 



(3.12) 



[0061] Eine GroBe, die beim gewahlten Ansatz einen wesentiichen EinfluB auf den maximal moglichen Massenslrom 
mzy, besitzt, ist der Umgebungsdruck P^y. Aus diesem Grund kann nicht von einem konstanten Wert dieser GroBe 
ausgegangen werden, sondern.es erfolgt eine 'Anpassung in der nachfolgend beschriebenen Art und:Weise. 
[0062] Der Wert des Umgebungsdruckes P J wird verandert, jj/enn der Betrag der KorrekturgroBe AApf^Q eine be- 
stimmte Schwelle uberschreitet oder wenn das Druckverhaltnis-^ groBer als eine wahlbare Konstante ist. Damit wird 
gewahrleistet, daB sowohl im Teil- als auch im Votlastbereich ein^ Umgebungsdruckanpassung erfolgen kann. 
[0083] Im folgenden wire! ein Mpdellabgleich fur luftmassengef uhrte Motorsteuerungssysteme erklart. Fur dieses 
System kann die in Figur 3 dargestellte Modellstruktur angegebep werden. 

[0064] Der D'rosselklappenstellungsfuhler 14 (Figur 1) liefert ein dem Offnungsgrad der Drosselklappe 11 entspre- 
chendes Signal, z.B. einen Drosselklappenoffnungswinkel. In einem Kennfeld der elektronischen Motorsteuerungs- 
einrichtung sind zu verschiedenen Werten dieses Drosselklappenoffnungswinkels zugehorige Werte fOr den reduzier- 
ten Querschnitt der Drosselklappe A^^p abgespeichert. Diese Zuordnung wird durch den Block "statisches Model!" in 
Figur 3 und in Figur 4 reprasentiert. Das Teilsystem "Saugrohrmodell" in den Figuren 3 und 4 reprasentiert das durch 
(2.7) beschriebene Verhalten. FOhrungsgroBe dieses Modellregelkreises ist der MeBwert des uber ein Segment ge- 
mittelten Luftmassenstromes an der Drosselklappe mm_LAM/. Wird als Regler in diesem Modellregelkreis ein Pl-Regler 
eingesetzt, so ist die bleibende Regelabweichung null, d.h. ModellgroBe und MeBgroBe des Luftmassenstromes an 
der Drosselklappe sind identisch. Die Pulsationserscheinungen des Luftmassenstromes an der Drosselklappe, die vor 
allem bei 4-Zylindermotoren zu beobachten sind, fuhren bei betragsbildenden Luftmassenmessern zu erheblichen 
positiven MeBfehlern und somit zu einer stark tehlerbehafteten FOhrungsgroBe. Durch eine Abschaltung des Reglers, 
d.h. einer Verkleinerung der Reglerparameter kann zum gesteuerten modellgestOtzten Betrieb ubergegangen werden. 
Bereiche, in denen die genannten Pulsationeri auftreten, konnen somit mit dem selben Verfahren unter Berucksichti- 
gung dynamischer Zusammenhange behandelt werden^ wie diejenigen Bereiche, in denen eine nahezu ungestorte 
FOhrungsgroBe vorliegt. Im Gegensatz zu Verfahren, die relevante MeBwerte nur In stationaren Betriebspunkten, be- 
rOeksichtigen, bleibt das beschriebene System nahezu uneingeschrankt arbeitsfahig. Bei Ausfall des Luftmassensi- 
gnals oder des Signals des Drosselklappensteliungsf Ohlers ist das vorgestellte System in der Lage, ein entsprechendes 
Ersatzsignal zu bilden; Bei Ausfall der FOhrungsgroBe muB der gesteuerte Betrieb realisiert werden, wahrend im an- 
deren Fall der geregelte Betrieb die kaum beeintrachtigte Funktionsfahigkeit des Systems garantiert. 
[0065] Der Block "Saugrohrmodell" reprasentiert die Verhaltnlsse wie sie anhand der Gleichung (2.7) beschrieben 
sind und hat demzufolge als AusgangsgroBe die ModellgroBe P^ sowie die zeitliche Ableilung^, und die GroBe 
Nach der Modellierung des SehsorObertragungsverhaltens d.h. des Uberlragungsverhaltens des Luftmassenmessers 
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und der Abtastung wird die ModellgroRe 'Lo^_u^ einer Mittelung unterzogen. so daB die gemittelte GroBe i« und der 
vom Luftmassenmesser gemessene durchschnittliche Luftmassenstrom ^m u,,. einem Vergleicher zugefuhrt werden 
konnen. Die Differenz beider Signale bewirkt eine Anderung AA^^^q reduzierten Stromungsquer-schnittes A^^^q. 
so daB stationar und instationar ein Mo-dellabgleich erfolgen kann. 

[0066] Fur saugrohrdruckgefuhrte Motorsteuerungssysteme wird die in Figur 4 dargestellte Modellstruktur angege- 
ben, wobei gleiche Blocke wie in Figur 3 gleiche Bezeichnungen tragen. Ebenso wie bei dem luftmassengefuhrten 
Motorsteuerungssystem reprasentiert das Teilsystem "SaugrohrmodelT, das durch die Differenzengleichung (2.7) be- 
schriebene Verhalten. FuhrungsgroBe dieses Modellregelkreises ist der MeBwert des uber ein Segment gemittelten 
Saugrohrdruckes 5. Wird ebenfalls wie in Figur 3 ein Pl-Regler eingesetzt, so ist im stationaren Fall der MeBwert 
des Druckes im Saugrohr P^^^mit der Modellgr6Be>,., idenlisch. Wie oben beschrieben, bleibt auch das vorliegende 
System nahezu uneingeschrankt arbeitsfahig, da bei Ausfall des Saugrohrdrucksignales oder des MeBwertes fur den 
Drosselklappenwinkel ein entsprechendes Ersatzsignal gebildet werden kann. 

[0067] Die durch das Saugrohrmodell erhaltenen ModellgroBen Pq, werden einem Block "Pradiktion" zugefuhrt. 
Da mit den Modellen auch die Druckanderungen im Saugrohr berechnet werden, konnen diese Druckanderungen dazu 
verwendet werden; den zukunftigen Druckverlauf im Saugrohr und damit die Zylinderluftmasse fur das nachste [A/+1] 
oderfur die nachsten Segmente [A/+H] zu schatzen. Die GroBe m2y/bzw. die GroBe m2y/[A/+1 ] dienen dann zur exakten 
Berechnung der Einspritzzeit, wahrend derer Kraftstoff eingespritzt wird. 



Patentanspruche 

1. Verfahren zum Bestimmen der einstromenden Luftmasse in den bzw. die Zylinder einer Brennkraftmaschine mit 

einem Ansaugsystem, das ein Saugrohr (10) und eine darin angeordnete Drosselklappe (11). sowie einen 
den Offnungsgrad der Drosselklappe (19) erfassenden Drosselklappenstellungsfuhler (14) aufweist, 
einem, ein Lastsignal ( ; P5 5) der Brennkraftmaschine erzeugenden Sensor (12; 13) 
einer elektrischen Steuerungseinrichtung, die auf der Grundtage des gemessenen Lastsignals ( 'mm_Lxm \Pss) 
und der Drehzahl der Brennkraftmaschine eine Grundeinspritzzeit berechnet, wobei 

die Verhaltnisse im Ansaugsystem mittels eines Saugrohrfullungsmodells nachgebildet werden, wobei als 
EingangsgroBen des Modells der Offnungsgrad der Drosselklappe (11), der Umgebungsdruck {Py) und die 
Ventilstellung reprasentierende Parameter herangezogen werden, 

eine ModellgroBe fur den Luftmassenstrom (^oa) an der DrosselWappe (11) mit Hilfe der DurchfluBgleichung 

idealer Gase durch Drosselstellen beschrieben wird, 

eine ModellgroBe^fur den Luftmassenstrom (mzw) in den bzw. in die Zylinder (17) als lineare Funktion des 

Saugrohrdrucks (P^) durch eine fvlassenbilanz der Luftmassenstrome (^oa, ^2^) beschrieben wird 

diese ModellgroBen uber eine Differentialgleichung verknupft werden, daraus^als bestimmende GroBe zur 

Ermittlung der tatsachlichen Last der Brennkraftmaschine der Saugrohrdruck (Pg) berechnet wird und 

aus dem linearen Zusammenhang zwischen berechnetem Saugrohrdruck (Pg) und der ModellgroBe fur den 

Luftmassenstrom (mz.*) in den bzw. in die Zylinder (17) durch Integration die in den bzw. in die Zylinder (17) 

einstromende Luftmasse {mzyi) erhalten wird. 

2. Verfahren nach Anspruch 1 = 
dadurch gekennzeichnet, 

daB das vom Lastsensor (12; 13) gemessene Lastsignal ( i.o*.t«M ; Ps_s) Korrektur und damit zum Abgleich der 
ModellgroBen (mz^J^in einem geschlossenen Rgelkreis herangezogen wird, wobei das Lastsignal ( : Ps_s) 
als FuhrungsgroBe des Regeikreises dient. 

3. Verfahren nach Anspruch 2, 
dadurch gekennzeichnet, 

daB der Abgleich im stationaren und/oder instationaren Betrieb der Brennkraftmaschine durchgefuhrt wird und 
dabei das Ubertragungsverhalten des Lastsensors (12; 13) berucksichtigt wird. 

4. Verfahren nach Anspruch 2, 
dadurch gekennzeichnet, 

daB jedem gemessenen Wert des Drosselklappenoffnungsgrades ein Wert eines reuzierlen Querschnitts der Dros- 
selklappe (Api^q) zugeordnet ist und der Abgleich der ModellgroBen durch Korrektur des reduzierten Querschnitts 
(^red) durch eine KorrekturgroBe (i^Api^Q) derart erfolgt, daB die Regelabweichung zwischen FuhrungsgroBe und 
entsprechender ModellgroBe minimiert wird. 
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5. Verfahren nach Anspruch 4, 
dadurch gekennzeichnet, 

daB der reduzierte Querschnitt (A^^^) stationaren Messungen am MotorprOfstand ermittelt wird und in einem 
Kennfeld eines Speichers der elektrischen Steuerungseinrichlung abgelegt ist. 
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6. Verfahren nach Anspruch 1, 
dadurch gekennzeichnet, 

daQ bei der Darstellung der ModellgroBe fur den Luftnnassenstrom (moK-) an der Drosselklappe (11) eine in der 
Durchflu3gleichung (Glelchung 2.2) vorhandene DurchfluBfunktlon (\^) in einzelne Abschnitte (i = 1..,k) unterteilt 
wird und diese Abschnitte durch Geradenabschnltte angenahert warden, wobei fur die Steigung (m|) und^ur das 
Absolutglied (nj) der jeweiligen Geradenabschnltte als Funktion des Verhaltnisses von Saugrohrdruck (Pq) und 
Umgebungsdruck (P^/) bestimmt sind und in einem Kennfeld abgelegt sind. 
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7. Verfahren nach Anspruch 1, 
dadurch gekennzeichnet, 

da3 die Steigung {y{) uiid das Absolutglied (yq) der linearen Funktion fur die ModellgroBe fur den Luftmassenstrom 
in den bzw. ind die Zylinder (m^w ) abhangig von mindestens einem der Parameter, Drehzahl der Brennkraftmaschi- 
ne, Zylinderzahl,.Saugrohrgeometrie, Temperatur der Luft (Tg) im Ansaugrohr (10) und Ventilsteuerzeichen fest- 
gelegt sind. 

8. Verfahren nach Anspruch 7, 
dadurch gekennzeichnet, 

daB die Parameter duch stationare Messungen am MotorprOfstand ermittelt werden und in Kennfelder abgelegt 
sind. 

9. Verfahren nach Anspruch 1,. 
dadurch gekennzeichnet, 

daB die in den Zylinder einstromende Luftmasse (/TJ2y/) durch die Beziehung 
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mzyi [ ~ 1] -h mzyi [ N\ 



berechnet wird, mit 

Ty^: Abtastzeit oder Segmentzeit 
^0 mz„ [A/1 ModellgroBe des Luftmassenstromes wahrend des aktuellen Abtastschrittes oder Segments 

«2v, [A/- 1] ModellgroBe des Luftmassenstromes wahrend des vergangeneri Abtastschrittes oder Segmentes. 
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10. Verfahren nach Anspruch 1 
dadurch gekennzeichnet, 

daB die in den bzw. in die Zylinder einstromende Luftmasse (m^y/) fur einen, bezuglich der aktuellen Lasterfassung 
zum Abtastzeitpunkt.(A/| in der Zukunft liegenden, bestimmten Pradiktionshorizont (H) geschatzt wird durch Schat- 
zung des entsprechenden Druckwertes nach folgender Beziehung: 
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mit 



T^: Abtastzeit oder Segmentzeit 

H: Pradiktionshorizont. Anzahl der in der Zukunft liegenden Abtastschritte 



1 \ 
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Yi : Steigung der llnearen Gleichung 
Yo- Absolutglied zur Beslimmung von m^y/ 
N: aktueller Abtastschritt 

11. Verfahren nach Anspruch 10, 
dadurch gekennzeichnet, 

da3 die Anzahl (H) von Segmenten, fur die das Lastsignal fur die Zukunft geschatzt werden soil, drehzahlabhangig 
festgelegt ist. 

Claims 

1. Method for determining the air mass flowing into the cylinder or cylinders of an internal combustion engine with 

'5 - an induction system which comprises an induction manifold (10) and a throttle valve (11 ) arranged within it as 

well as a throttle valve position sensor (14) which detects the degree of opening of the throttle valve (19), 

a sensor (12; 13) which produces a load signal ( iDA^^wv ; Ps_s^ internal combustion engine, 

20 - an electrical control device, which calculates a basic injection time on the basis of the measured load signal 

( : Ps_s) speed of the internal combustion engine, whereby 

the conditions in the induction system are simulated by an induction manifold charging model, wherein pa- 
rameters representing the degree of opening of the throttle valve (11 ), the ambient pressure (P^y) and the valve 
2S position are used as input variables of the model, 

a model variable for the air mass flow {^^k ) at the throttle valve (11 ) is described by means of the flow-through 
equation of ideal gases through throttling points, 

30 - a model variable for the air mass flow {'mz.i) into the cylinder or cylinders (17) is described as a linear function 

of the induction manifold pressure (Pq) by a mass balance of the air mass flows (moK, ^z>0. 

these model variables are linked by a differential equation, the induction manifold pressure {Pq) being calcu- 
lated therefrom as a defining variable for determining the actual load of the internal combustion engine, and 

35 

the air mass (mzrt) flowing into the cylinder or cylinders (17) is obtained from the linear relationship between 
calculated induction manifold pressure (P^) and the model variable for the air mass flow (m^^/) in the cylinder 
or cylinders (17) by integration. 

40 2. Method in accordance with claim 1, characterized in that the load signal {'mnK,^^^: Ps_s) measured by the load 
sensor (12; 13) is used to correct and therefore to adjust the model variables (mw) in a closed control circuit, the 
load signal ( 

i Ps s) serving as reference variable of the control circuit. 

3. Method in accordance with claim 2, characterized in that the adjustment is made in the steady-state and/or transient 
45 mode of the internal combustion engine, allowing for the transmission behaviour of the load sensor (12; 13). 

4. Method in accordance with claim 2, characterized in that a value of a reduced cross-section of the throttle valve 
('^red) allocated to each measured variable of the degree of opening of the throttle valve and the model variables 
are adjusted by correcting the reduced cross-section (>Ar£q) by a correction variable (A^Ar^d) such that the deviation 

so between reference variable and corresponding model variable is minimized. 

5. Method in accordance with claim 4, characterized in that the reduced cross-section (^Apj^p) is determined from 
steady-state measurements on the engine test stand and is entered into a map in a memory in the electric control 
device. 
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6. Method in accordance vi/ilh claim 1 , characterized in that, when portraying the model variable for the air mass flow 
(rno,. ) at the throttle valve (11). a flow-through function (\}/) present in the flow-through equation (equation 2.2) is 
divided into individual sections (i = 1...k) and these sections are approximated by straight line sections, the re- 
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.spective straight line sections being determined as a function of the ratio of induction manifold pressure {Pq) to 
ambient pressure {PJ) for the gradient (mj) and for the absolute term (n|) and being entered into a map. 

7. Method in accordance with claim 1 . characterized in that the gradient (Yi) and the absolute term (jq) of the linear 
function for the model variable for the air mass flow into the cylinder or cylinders i^zyi ) are established as a function 
of at least one of the parameters of speed of the internal combustion engine, number of cylinders, induction manifold 
geometry, temperature of the air (Tg) in the induction manifold (10) and valve control symbol. 

8. Method in accordance with claimi 7, characterized in that the pararneters are determined by steady -state meas- 
urements on the engine test stand and are entered into maps. 

9. Method in accordance with claim 1 , characterized in that the air mass (n?2y/) flowing into the cylinders is calculated 
by the equation 
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where 

Izy^ [Ad 

I^.[A/-1] 



is the sampling interval or segment time 

is the model. variable of the air mass flow during the current sampling step or segment 
is the model variable of the air mass flow during the previous sampling step or segment. 



10. Method in accordance with claim 1 , characterized in that the air mass (m^yj) flowing into the cylinder or cylinders 
is estimated for a defined prediction horizon (H) located In the future with respect to the current load detection at 
the moment of sampling [A/] by estimating the corresponding pressure value in accordance with the following 
equation: 
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Ps W + (if - 0. 5) . ^ . [P, [i^ - 1] + P, [ A^]l 
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where 

H: 
Yi : 

Vo- 

N: 



sampling interval or segment time 

prediction horizon, nunnber of sampling steps located in the future 
gradient of the linear equation 
absolute term for determining 
current sampling step 



11. Method in accordance with claim 10, characterized in that the number (H) of segments for which the load signal 
is to be estimated for the future is established as a function of the speed. 



55 Revendications 

1. Procede de determination du debit massique d'air inlroduit dans le ou les cylindres d*un moteur a combuslion 
interne, comprenant : 
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un systeme d'admission comportant un collecteur d'admission (10) et un papillon des gaz (11) dispose dans 
celui-ci, ainsi qu'un dispositif capteur de position de papillon des gaz (14) detectant le degre d'ouverture du 
papillon des gaz (19), 

un capteur (12 ; 13) produisant un signal de charge ( lot.utu ; P5 s) moteur a combustion interne, 

un dispositif de commande electronique qui calcule une duree d'injection de base a partir du signal de charge 

mesure ( iw.ut>« ; Ps s) ®^ ^® vltesse de rotation du moteur a combustion interne, 

selon lequel : 

les conditions dans le systeme d'admission sont simulees au moyen d'un modele de remplissage de collecteur 
d'admission, le degre d'ouverture du papillon des gaz (11), la pression ambiante {PJ) et des parametres re- 
presentant la position de soupapes etant utilises en tant que grandeurs d'entree du modele, 
une grandeur de modele correspondant au debit massique d'air (^dj:) a I'endroit du papillon des gaz (11) est 
decrite au moyen de I'equation d'ecoulement des gaz parfalts a travers des emplacements d'etranglement, 
une grandeur de modele correspondant au debit massique d'air (mzw) dans le ou les cylindres (17) est decrite 
comme fonction lineaire de la pression de collecteur d'admission {Pq) au moyen d'un bilan des debits massi- 
ques d'air {Iok, -2.0, 

ces grandeurs de modele sont liees au moyen d'une equation differentieile et, a partir de celle-ci, la pression 
de collecteur d'admission (PJ est calculee comme grandeur decisive pour la determination de la charge ef- 
fective du moteur a combustion interne et 

A 

a partir de la relation lineaire se presentant entre la pression de collecteur d'admission calculee (Pg) et de la 
grandeur de modele correspondant au debit massique d'air («^w) dans le ou les cylindres (17), le debit d'air 
(^Zy/) introduit dans le ou les cylindres (17) est obtenu par integration. 

Procede suivant la revendication 1 , caracterise en ce que le signal de charge ( moi.-.i.wu ; P^ 5) mesure par le capteur 
de charge (12 ; 13) est utilise dans un circuit de regulation ferme pour la correction et done pour la compensation 
des grandeurs de modele («?« ), le signal de charge ( ^„r_t.«.« ; Ps_s) servant de grandeur de pilotage du circuit de 
regulation. 

Procede suivant la revendication 2. caracter-ise en ce que la compensation est effectuee en regime permanent et/ 
ou en regime non permanent du moteur a combustion interne et qu'il est de ce fait tenu compte du comportement 
de transmission du capteur de charge (12 ; 13). 

Procede suivant la revendication 2, caracterise en ce qu'une valeur d'une section transversale reduite (>A^5d) du 
papillon des gaz est associee a chaque valeur mesuree du degre d'ouverture de papillon des gaz et que la com- 
pensation des grandeurs de modele par correction de la section transversale reduite (Ap^p) au moyen d'une 
grandeur de correction [AAf^^p) est effectuee de fagon telle que Tecart de regulation entre la grandeur de pilotage 
et la grandeur de modele correspondante est rendu minimal. 

Procede suivant la revendication 4, caracterise en ce que la section transversale reduite (^A^^q) est determinee 
a partir de mesures effectuees en regime permanent sur banc d'essais de moteur et est rangee en memoire dans 
une table caracteristique d'une memoire du dispositif de commande electrique. 

Procede suivant la revendication 1. caracterise en ce que, pour la representation de la grandeur de modele cor- 
respondant au debit massique d'air (^o*. ) a I'endroit du papillon des gaz (11 ), une fonction d'ecoulement {x^f) existant 
dans I'equation d'ecoulement (equation 2.2) est divisee en segments individuels (i = 1 . . k) et ces segments font 
I'objet d'une approximation au moyen de segments de droite, la pente (mj) et I'abscisse a I'origine (n;) des^segments 
de droite respectifs etant determines en fonction du rapport de la pression de collecteur d'admission (P^) et de la 
pression ambiante (Py) et etant rangees en memoire dans une table caracteristique. 

Procede suivant la revendication 1, caracterise en ce que la pente (y^) et I'abscisse a I'origine (yq) de la fonction 
lineaire correspondant a la grandeur de modele prevue pour ie debit massique d'air {m^yi) dans le ou les cylindres 
sont determinees en fonction d'au moins I'un des parametres que sont la vitesse de rotation du moteur a combustion 
interne, le nombre de cylindres. les parametres geometriques du collecteur d'admission, la temperature (T^) de 
I'air dans le collecteur d'admission (10) et les durees de commande de soupapes. 

Procede suivant la revendication 7. caracterise en ce que les parametres sont determines au moyen de mesures 
effectuees en regime permanent sur banc d'essais de moteur et sont ranges en memoire dans des tables carac- 
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teristiques. 



Procede suivant la revendication 1 , caracterise en ce que le debit masslque d*air {m^yl^ penetrant dans le cylindre 
est calcule au moyen de la relation, 



10 



IS 



f ^ 



mzy{[N - 1] + mzyi [N] 



designant la duree d'analyse par echantillonnage ou segment de temps, 

la grandeur de modele du debit massique d'air pendant le pas d'analyse ou segment actuel, 

la grandeur de modele du debit massique d'air pendant le pas d'analyse ou segment passe. 



10. Procede suivant la revendication 1, caracterise en ce que le debit masslque d'air {m^yi) penetrant dans le ou les 
cylindres est estime pour un horizon de prediction (H) determine, situe dans le futur par rapport a la determination 
de charge actuelle a I'instant d'analyse par echantillonnage [A/], par estimation de la valeur de pression corres- 
pondarite conformement a la relation suivante : 



20 



25 



avec 



30 



35 



H 

Ti 
Yo 
N 



duree d'analyse par echantillonnage ou segment de temps, 

horizon de prediction, nombre des pas. d'analyse par echantillonnage se trouvant dans le futur, 
pente de I'equation lineaire, 

abscisse a Torigine servant a la determination de mj;^, 
pas d'analyse par echantillonnage actuel. 



11. Procede suivant la revendication 10, caracterise en ce que le nombre (H) de segments pour lesquels le signal de 
charge doit etre estime dans le futur est determine en fonction de la Vitesse de rotation. 



40 



45 



so 
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